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Abstract 

The long range properties of four-dimensional compact U(l) lattice gauge theory 
with the Wilson action in the confinement phase is studied by using the multi- 
level algorithm. The static potential, force and flux-tube profile between two static 
charges are successfully measured from the correlation function involving the Polyakov 
loop. The universality of the coefficient of the 1/r correction to the static potential, 
known as the Liischer term, and the transversal width of the flux-tube profile as a 
function of its length are investigated. While the result supports the presence of the 
1/r correction, the width of the flux tube shows an almost constant behavior at a 
large distance. 



1 Introduction 



There is a conjecture that (non supersymmetric) confining gauge theories in 
the infrared regime are related to an effective bosonic string theory. If this 
happens, the asymptotic behavior of the potential between static charges sep- 
arated by distance r is expected to be parametrized as 

V{r)=(Tr + fi + ^ + 0{\). (1) 



Here a is the string tension, which characterizes the strength of the confining 
force of static charges, and /i denotes a constant. The third term is the zero 
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point Casimir energy for an open bosonic string with fixed boundary. This 
correction is known as the Liischer term [1] and the coefficient 7 is considered 
to be universal, such that it does not depend on the gauge group but only on 
the space-time dimension d through 7 = —7r{d — 2)/24. The effective bosonic 
string theory also predicts that the width of the field energy distribution of 
the flux tube diverges logarithmically as r — > 00 [2]. Recent Monte Carlo 
simulations of various lattice gauge theories support the universality of 7 in 
Eq. (1) with high accuracy: the confinement phase of lattice gauge theory 
(LGT) in 3D [3,4,5], SU(2) LGT in 3D [6,7] and in 4D [8], SU(3) LGT in 3D [9] 
and in 4D [9,10]. Moreover, in Refs. [6,10], the spectrum of the string states 
have been computed, which further support the effective string description of 
confining gauge theories. 

In this context, we are now interested in the 4D compact U(l) LGT with 
the Wilson action. This theory possesses a confinement phase analogous to 
non-Abelian gauge theories below the critical coupling /3 < /3c ~ 1.011 (pre- 
cise value can be found in Ref. [11]). This is due to the presence of magnetic 
monopoles [12], which cause the dual Meissner effect like in a dual supercon- 
ductor [13,14] when electric sources are introduced in the vacuum; the electric 
fiux is squeezed into a fiux tube by the circulating monopole supercurrent, 
which leads to a linear rising potential between static charges. In fact, the 
measurements of the U(l) flux-tube proflle in the conflnement phase have 
been reported in Refs. [15,16], which support the dual superconductor pic- 
ture. To answer the questions i) whether the static potential in this theory 
also contain the universal correction, ii) how the width of the fiux-tube profile 
behaves as a function of r, in this paper, we investigate the long range proper- 
ties of the potential, force and the flux-tube proflle between two static charges. 
Here we are going to use the Polyakov loop correlation function (PLCF: a pair 
of Polyakov loops separated by a distance r) as an external source. Contrary 
to the use of the Wilson loop {W{r, t)), if we use the PLCF {P*P{r)), we do 
not need to care about the t dependence of the results and can extract the the 
ground state easily as long as the lattice temporal extension is large enough. 
This is important because most of analytical predictions are given for such 
a ground state. However, the measurement of the PLCF in the conflnement 
phase with large r is a quite difficult task since the expectation value becomes 
exponentially smaller with increasing r. Moreover, due to the strong coupling 
nature of the theory, the signal-to-noise ratio is very small from the begin- 
ning. In principle, one would need enormous statistics and computation time 
to identify such small expectation values. 

Recently, Liischer and Weisz (LW) have proposed a multi-level algorithm for 
pure LGT [17] to compute the expectation value of a Wilson loop for a large 
size and a PLCF for large r with exponentially reduced statistical errors. 
They noted that its algorithmic idea is essentially the same as in the multi-hit 
method [18] but is applied to pairs of links instead of single finks. In fact, based 
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on their algorithm they have confirmed the presence of an universal potential 
7/r in SU(3) LGT [9]. The LW algorithm is applicable to other LGTs as long 
as a local gauge action is simulated. The Wilson action is the easiest gauge 
action to adopt. The studies of 3D SU(2) LGT [6,7] also use this algorithm. 
We then expect that the LW algorithm can also be apphed to compact U(l) 
LGT, which will help to overcome its numerical difficulties. In the original 
work [9], although the potential and its derivatives with respect to r, force 
etc., have been of interest, we find that this is also applicable to measuring 
the flux-tube proflle as well as the glueball mass measurements [19]. 

In section 2, we describe the LW algorithm for the measurements of the static 
potential and force from the PLGF. We also explain how to measure the flux- 
tube proflle in this context. In sections 3 and 4, we present simulation details 
and the numerical results on the PLCF/potential/force and on the the flux- 
tube proflles, respectively, where several analyses of the data are given. The 
section 5 is a summary. A part of these studies has been presented at the 
Lattice 2003 at Tsukuba, Japan [20]. 



2 Numerical procedures 

In this section we describe how to measure the static potential, force and 
flux-tube proflle with the LW algorithm. 

The Wilson gauge action of the compact U(l) LGT is given by 

S[U]^f3^^{l-Re[U,.{m)]} , (2) 

where U^viiT^) G U{1) are plaqucttc variables constructed from link variables 
Un{m) = exp(i6'^(m)) e U{1). We consider the lattice volume A''^ x A^^^, where 
Nt is an even number. The four-dimensional sites are labeled by m = (m^, mj) 
with 777,5 — (?77i, 7772, 7773). Periodic boundary conditions are imposed for all 
directions. A sequence of independent gauge fleld conflgurations, labeled by 
ic = 1,2, Nc, are generated by using a mixture of heatbath (HB) and 
over-relaxation (OR) link updates as usual. 

2.1 The static potential and force from the PLCF 

We adopt the terminology of Ref. [17]. To measure the PLGF {P*P{R)) for 
a charge distance R — r / a, we flrst take sub-lattice averages of the two-hnk 
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l=Nt+I 

Jf"T] 

[p*p] = — !— y • 



m m + Ri 

Fig. 1. How to construct the [-P*-P]ic with the LW algorithm. [• • •] denotes the 
sub-lattice average. 

correlators 

T{m;R;i) = Ul{m)U4{m + Ri) , (3) 

as 

T^'^\ms,fnt;R;i) = [T{m; R;i)T{m + A; R;i)] , (4) 

where i — 1,2,3 are possible directions of two static charges and fht — 
1,3, Nt — 1. The sub-lattice average is achieved by updating link vari- 
ables (with a mixture of HB/OR) except for the spatial links at the time slice 
fht- We call this procedure the internal update. We repeat the internal update 
until reasonably stable values for T^^^ are obtained. 

Then, the PLCF at a spatial site is constructed from T*^^^ as 

P*P{ms; R; i) = P*{ms)P{ms + Rii) 
= Re T(2) (m„ 1; R; i)!^^) (m^, 3; i?; i) • • • T^^^ (m^, Nt - 1; R; i) . (5) 

We take the average with respect to and i, which provides the value of the 
PLCF for icth configuration, [P*P(i?)]-^. For a schematic understanding, see 
Fig. 1. The desired expectation value {P*P{R)) is calculated from the average 
of [P*P(i?)].^ for z, = l,2,...,7V,. 

The static potential and the corresponding force are taken as (neglecting terms 
of 0(e-(^^)^*)) 

aV{R) = -^HP*P{R)) , (6) 
a^F{R) = aV{R) - aV{R - 1) , (7) 
where R^R- 1/2. 
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In the actual measurements, we have also applied the multi-hit technique [18] 
to the timelike link variables Ui(rn) for i? > 2 before constructing the two- 
link correlators (3). It was further helpful to reduce the statistical errors of 
the PLCF within a hmited CPU time. In U(l) LGT, this procedure is given 
by the replacement 

^ imW,{m)\) W,{m) 



where Jq and /i denote the zeroth- and first-order modified Bessel functions 
and Wfj,{m) the sum of six staples around U^{m): 



Some comments associated with the LW algorithm are in order. The number 
of the internal updates Mupd is an optimization parameter that has to be 
tuned for efficient performance of the computation. If one wants to compute 
the force, which requires two values of the potential at different r, it is useful to 
compute all i? = 1,2,..., i?max in one run without changing N^^-p^ depending 
on i?, although a small number of A^iupd is enough for a short distance. This 
is because data among different i?'s are highly correlated, which leads to a 
significant cancellation of the statistical errors in the difference [9] . Practically, 
one may regard not only the PLCF but also the potential and force as the 
primary observables and apply the jackknife analysis for the evaluation of the 
statistical errors. 



2.2 The flux-tuhe profile 



In order to measure the flux-tube proflle, one needs to compute a correlation 
function of the type 



iO(n)), ^ - . (10, 

where 0{n) is a local operator, (• ■ ■)j denotes an average in the vacuum with 
the PLCF, and (• • ^ an average in the vacuum without such a source. For a 
parity-odd local operator, we do not need the second term since it gives no 
contribution, (O)o = 0. However, if one is interested in a parity-even local 
operator such as the action density cos 6*^,^(71), where 9n^{n) is the phase of 
U^i,{n), one needs to subtract out the vacuum expectation value. It is noted 
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l=Nt+l 



Nt-1 



[P*PO] = 



3Ns3(Ni/2) 



m n m + Ri 



Fig. 2. How to construct \P*PO\i^ with the LW algorithm. [• • •] denotes the 
sub-lattice average and the square represents a local operator. 

that to receive maximum benefit from the LW algorithm, the parity-odd local 
operator is preferable [19]. 

To measure (P*P(9(n))o on the mid-plane between two static charges, we 
parameterize the position of the local operator n as 



n — m^ {R/2)i + xj + yk , 



where i is the direction of two static charges and j — k specify a 2D plane 
perpendicular to i. By constructing the two- link- local-operator correlators 



0(m; n; R; i) = Ul{m)Ui{m + Ri)0{n) , 

we compute sub-lattice averages of the correlation function 
TO^^) (m^, fht; n; R; i) = [T(m; R] i)<0{m + 4; n; R; i)] . 



;i2) 



(13) 

Combining TO*^^^ and T^^^ in Eq. (4), we obtain the PLCF involving a local 
operator at site as 

P*PO{m,;x,y;R;i) 

^ Re |TO(2)(m5, 1; n; R; i)T^^\ms, 3;R;i)--- T^^) (m„ Nt - 1; R; i) + 

rUs^Nt - l;n;R;i) \ . (14) 



iNt/2) 

h T^^^ {rus, 1; R; i)T^'^^ {rus, 3; R; i) 



The average with respect to rris and i provides [P*PO{x^ y; R)] - . For a schematic 
understanding, see Fig. 2. Wc repeat the same procedure as for the PLCF to 
get the final expectation value {P*PO{x,y; R)) . 

As a local operator, we use the electric field operator (parity odd) as 

OE{n) = i9^u{n) = i{9^y{n) - 27rn/,^(n)) , (15) 
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where O^uin) G [—4:71, An] and n^y{n) G [0, ±1, ±2] is the modulo 27r of 6^y{n), 
which corresponds to the magnetic Dirac string. Hence, one has 9^y{n) G 
[— TT, tt]. Moreover, we use the monopole current operator (parity odd) to detect 
the circulating monopole supercurrent as 



Here is defined as the boundary of the magnetic Dirac string [12] as 



n denotes the dual site n + (i + 2 + 3 + 4)/2. 

We have chosen these local operators because of the possibihty to relate the 
U(l) flux tube and the classical flux-tube solution of the dual Ginzburg- 
Landau (DGL) theory. In the DGL theory, the circulating monopole supercur- 
rent induces the solenoidal electric field through the dual Ampere law, which 
plays a role in cancelling the Coulombic field induced by static charges at a 
distant place. In this sense the measurement of the monopole current profile 
is useful to judge whether the total electric flux is indeed squeezed or not. In 
Appendix A, we show a numerical evidence of such a cancellation mechanism 
of the electric flux inside the U(l) flux tube. We may call this the composite 
structure of the U(l) flux tube. 

Note that the definitions of the electric field and monopole current operators 
are not unique. For instance, Cheluvaraja et al. [21], have proposed an alterna- 
tive definition to satisfy the the Maxwell equations at finite lattice spacing. It 
would also be interesting to study how the measured flux-tube profiles change 
with their operators. 



3 Numerical results : Static potential and force from the PLCF 

In this section, we first present the simulation details associated with the LW 
algorithm for the measurements of the potential and force from the PLCF, 

and then, we show the corresponding numerical results. Some analyses are 
also performed for the potential and force; the potential is fitted with several 
ansatze. The behavior of the force is compared with the function derived from 



(16) 



kn{n) = T.^nvapduriaiiin + fi) G [0, ±1, ±2] . 



(17) 



Eq. (1). 
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3.1 Simulation details 



We use a 16^ lattice. The (3 values, the number of internal updates Mupd, the 
number of configurations A^'conf, details of one Monte Carlo update (HB/OR), 
and the range of measured distance between static charges r/a, are summa- 
rized in Table 1. Although wc have not checked the finite volume effect, as 
we will see later, our lattice volume itself is reasonably large even near the 
phase transition point: (~ 3.5ro)^ at /3 = 1.01. In addition to this, we restrict 
ourselves to measure the potential up to r/a = 6. 

Table 1 

Parameter setting 



/3 


-^iupd 


-A^conf 


HB/OR 


r/a range 


0.98 


10000 


1050 


1/3 


1-6 


0.99 


8000 


1250 


1/3 


1-6 


1.00 


5000 


2000 


1/3 


1-6 


1.005 


3000 


2400 


1/5 


1-6 


1.01 


1000 


3200 


1/5 


1-6 



In Fig. 3, we show an example how wc have optimized A^mpd depending on (5. 
This plot shows a typical behavior of a PLCF at /? = 0.98 for one configuration 
[P*P{r/ay\ as a function of iViupd. This figure tells us that, for instance, if we 
are interested in up to r/a = 4, iViup^ = 1000 would be enough. However, if 
r/a = 6 is of interest, we may need to take Mupd > 8000. 

3.2 Static potential and force from the PLCF 

In Table 2, we summarize the expectation values of the PLCF, potential and 
force for all /3 values measured. Note that it is possible to identify signals of 
the PLCF even when (P*P) = 10~^ ~ 10~^^ with the la error varying from 
0.4 to 8 %. 

One may think that the investigation of the second derivative of the potential 
as in Refs. [6,9] for the direct identification of the coefficient of 1/r potential 
is also interesting. However we have not succeeded to get reliable data for 
this. The result was strongly dependent on the definition of the lattice second 
derivative. This may be due to the reason that in U(l) LGT the rotational 
invariance is not well-recovered compared to non-Abelian gauge theories even 
near the critical coupling. 
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Fig. 3. Typical behavior of [P*P] for various r at /3 = 0.98 as a function of A?iupd- 
When A^iupd is not sufficient, [P*P] often takes negative values during the internal 
update, typically for large r, where lines are broken. 
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Table 2 

The PLCF, static potential and force. 





r/a 


{P*P{ 


r/a)) 




aV{r /a) 


r/a 


a^F{f/a) 


0.98 


1 


1.799(3) 


X 10" 


-4 


0.5389(1) 


1.5 


0.4080(3) 




2 


2.627(14) 


X 10" 


-7 


0.9470(3) 


2.5 


0.3439(5) 




3 


1.072(14) 


X 10- 


-9 


1.2909(8) 


3.5 


0.3150(6) 




4 


6.96(15) 


X 10" 


1 o 
-LZ 


1.6058(14) 


4.5 


0.3024(12) 




5 


5.54(17) 


X 10" 


1 A 

-14 


1.9083(19) 


5.5 


0.3003(52) 




6 


4.77(35) 


X 10" 


-io 


2.2086(59) 






0.99 


1 


3.186(6) 


X 10" 


-4 


0.5032(1) 


1.5 


0.3562(3) 




2 


1.067(7) 


X 10" 


-6 


0.8594(4) 


2.5 


0.2872(4) 




3 


1.080(14) 


X 10" 


Q 

-o 


1.1466(8) 


3.5 


0.2585(7) 




4 


1.733(37) 


X 10" 


1 n 
-lU 


1.4051(14) 


4.5 


0.2474(12) 




5 


3.35(13) 


X 10" 


-12 


1.6525(24) 


5.5 


0.2444(40) 




6 


7.01(62) 


X 10" 


1 A 
-14 


1.8970(59) 






1.00 


1 


6.507(11) 


X 10" 


-4 


0.4586(1) 


1.5 


0.2920(2) 




2 


6.086(27) 


X 10" 


-6 


0.7506(3) 


2.5 


0.2193(5) 




3 


1.825(18) 


X 10" 


- i 


0.9699(6) 


3.5 


0.1907(14) 




4 


8.66(16) 


X 10" 


n 

y 


1.1606(12) 


4.5 


0.1797(27) 




5 


4.94(22) 


X 10" 


-lU 


1.3404(20) 


5.5 


0.1791(35) 




6 


3.00(23) 


X 10" 


-11 


1.5195(45) 






1.005 


1 


1.045(2) 


X 10" 


-3 


0.4290(1) 


1.5 


0.2504(2) 




2 


1.902(10) 


X 10" 


-5 


0.6794(3) 


2.5 


0.1770(4) 




3 


1.121(10) 


X 10" 


-D 


0.8565(7) 


3.5 


0.1498(6) 




4 


1.027(21) 


X 10" 


■7 
- / 


1.0063(13) 


4.5 


0.1389(12) 




5 


1.136(42) 


X 10" 


-8 


1.1452(22) 


5.5 


0.1350(35) 




6 


1.42(11) 


X 10" 


-9 


1.2802(52) 






1.01 


1 


2.152(8) 


X 10" 


-3 


0.3839(2) 


1.5 


0.1882(3) 




2 


1.061(10) 


X 10" 


-4 


0.5720(4) 


2.5 


0.1150(3) 




3 


1.687(24) 


X 10" 


-5 


0.6871(7) 


3.5 


0.0891(5) 




4 


4.076(79) 


X 10" 


-6 


0.7762(12) 


4.5 


0.0784(9) 




5 


1.180(37) 


X 10" 


-6 


0.8546(15) 


5.5 


0.0755(27) 




6 


3.79(22) 


X 10" 


-7 


0.9301(45) 
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3.3 Analysis 



In this subsection, in order to see the presence/ absence of the universal 7/r 
correction and of more higher order corrections to the static potential, we 
fit the static potential assuming the several explicit forms which are close to 
Eq. (1): 



C 

\/i(r) = ar + - , (18) 



r 



V^2(r) = (7r + At+^ , (19) 
r 

V,{r) = ar + ^^+^-{l+^^ , (20) 

where 7 = -7r/12 ~ 0.262. The form of is motivated by Ref. [9]. 

Before carrying out the fitting, the PLCF has been averaged in bins over inter- 
vals between 50 and 160 time units (representing 5000 and 16000 iterations) 
depending on {3 values to reduce the autocorrelation. The potential has been 
computed from the nested PLCF. For each fitting function, the means of the 
fitting parameters have been determined from the minimum of the which is 
defined with the covariance matrix so as to take into account the correlation 
among different r's. The errors of the fitting parameters have been estimated 
from the distribution of the jackknife samples of the fitting parameters. The 
fit range has been fixed so that the mean is consistent with that evaluated by 
using only the diagonal part of the covariance matrix. 

The fitting results are summarized in tables 3, 4 and 5. In Fig. 6 we have 
plotted the potential and the various fitting curves. As seen from this figure 
all curves are on the data and practically indistinguishable from each other 
apart from the point i? = 1, which lies outside of our fitting range. In all 
cases the orders of the x^/^df are one or smaller than one. The coefficient C 
obtained from Vi{r) is slightly different from the theoretical value of 7, but 
not by much. However this maybe due to the existence of the Coulombic 1/r 
potential in the current fitting region. An indication of this is seen in the fit 
by the potential V2(r) where we have had to drop the point R — 2. V^{r) gives 
a better fit, but of course it also contains an extra parameter. 

At this stage, one may wonder about the relevance of 1/r term at long dis- 
tances, since all investigated types of the function fit the potential very well. 
To test this we have fitted the potential with the simple form Vi{r) = or + 
with r = [3,6] and compared with the result of V2. It turned out that the 
fit was distinctly worse; the became 10 ~ 20 times larger than that of V2 
and the means were not consistent with that determined from the use of the 
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diagonal part of the covariance matrix. 

To compare tiie fit results among different p values, we introduce a scale 
To based on Sommer's relation F(ro)rQ = 1.65 [22]. We have found ro/a — 
2.13, 2.37, 2.83(1), 3.28(1), 4.60(3) for /3 = 0.98, 0.99, 1.00, 1.005 and 1.01, 
respectively. For first two f3 values there were no error at this order. Using 
tq, we plot oTq from Vi, V2 and V3 as a function of (3 in Fig. 4. The bottom 
axis corresponding to Vi and V3 are slightly shifted in the plot to distinguish 
them from each other. We see slight differences among the string tensions for 
different ansatze of the potential. The string tensions at /3 =0.98 and 0.99 
show a good scaling behavior with respect to Tq. For f3 > 1.00, o"jrQ start to 
grow, which suggests that the approach to the phase transition point of the 
string tensions and ro are different. In Fig. 5, we plot b/ro against P, where b 
is a coefficient of the 1/r^ potential in Eq. (20). It is interesting to find that 
at (3 =0.98 and 0.99, b/rQ seems to be saturated. As j3 increases, however, it 
falls down to zero. Since for large f3 we can look at only the short range of the 
potential, this result may suggest that b is irrelevant for such a range. 

We then show the potential for all /3 as a function of r/ro in Figs. 7. To subtract 
the constant /x in the potential plot, we have used the value obtained by the 
Vi fit. We find that the potential beautifully falls onto one curve except for 
the data from (3 = 1.01. The reason for this exception can be understood from 
Fig. 4, which shows that the string tension oTq grows as the (3 approaches to 
the phase transition point. We remark that the result was insensitive to the 
choice of the potential form used in the fit (corresponding figure from the V2 
fit is found in Ref. [20]). 

In Fig. 8, we plot the force for all as a function of r/ro and the expected 
function from the potential V2: -F2(r) = dV2{r)/dr = a — 7/r^. It should be 
noted that this function contains no fitting parameter, since Sommer's relation 
gives a fixed value for the string tension oTq = 1.65 — 7r/12 ~ 1.39. We find 
that the general behavior of the force data seems to be described by this 
function. Although there are shght differences between the curve and the data 
at long distance, we consider that this result supports the universality of 7/r 
correction to the static potential. To make a more precise statement, however, 
one has to control various systematic effects which enter in this analysis, for 
instance, the contribution of 1/r^ force which originates from the Coulombic 
electric field (partially discussed in Appendix A), and the possibility of higher 
order corrections to the static potential, which are not universal and vary from 
theory to theory. Surprisingly, another feature in common with non-Abelian 
gauge theories [23] is that this function also fits the data down to relatively 
short distance to r/ro ~ 0.3. For this, there is as yet no explanation. 
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Table 3 

Potential fit hy Vi{r) = ar + n + ^. 



p 


2 




C 


fit range (r/a) 


0.98 


0.286(1) 


0.546(5) 


-0.344(6) 


2-6 


0.99 


0.230(1) 


0.573(4) 


-0.346(4) 


2-6 


1.00 


0.162(1) 


0.598(3) 


-0.342(3) 


2-6 


1.005 


0.122(1) 


0.597(3) 


-0.326(3) 


2-6 


1.01 


0.0G19(.3) 


0.595(1) 


-0.305(1) 


2-G 



Table 4 

Potential fit by V2{r) = ar + fj, + ^ with 7 = — 7r/12. 



p 






fit range (r/a) 


0.98 


0.294(1) 


0.497(2) 


3-6 


0.99 


0.238(1) 


0.521(2) 


3-6 


1.00 


0.171(1) 


0.545(1) 


3-6 


1.005 


0.130(1) 


0.554(1) 


3-6 


1.01 


0.0709(8) 


0.565(1) 


3-6 



Tabic 5 

Potential fit by Vsir) = ar + fi + ^ {l + ^) with 7 = -7r/12. 



/5 




fia 


b/a 


fit range (r/a) 


0.98 


0.290(1) 


0.517(4) 


0.281(23) 


2-6 


0.99 


0.233(1) 


0.543(3) 


0.289(15) 


2-6 


1.00 


0.166(1) 


0.567(2) 


0.267(9) 


2-6 


1.005 


0.126(1) 


0.572(2) 


0.208(10) 


2-6 


1.01 


0.0667(5) 


0.578(1) 


0.138(5) 


2-6 
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4 Numerical results : Flux-tube profile 

In this section, we show numerical results on the flux-tube profile. We then 
investigate how the width of the flux tube behaves as a function of r based on 
the DGL analysis. 

4-1 Simulation details 

The P values, the lattice volume, and details of one Monte Carlo update 
(HB/OR) are the same as the measurement of the PLCF. Since in this case 
we do not compute derivatives with respect to r, we have changed Mupd de- 
pending on r in order to achieve a reasonable performance. The number of 
^iupd is summarized in Table 6. We have measured the profile on the mid- 
plane between charges as described in subsection 2.2. In order to compare the 
profile among different /3's and r's easily, we take the cylindrical average of 
the 2D profile; we define the radius p — -|- and the azimuthal angle 
around z axis as (/? = tan~'^(y/x). 

Table 6 

The number of A^iupd for the measurement of the flux-tube profile. We gave up the 
profile measurement for r/a = 6at/3 = 0.98 because of a practical reason. 



/3 


r/a 


-^iupd 


r/a 


-^iupd 


r/a 


^iupd 


r/a 


-^iupd 


0.98 


3 


200 


4 


1000 


5 


8000 


6 




0.99 


3 


200 


4 


1000 


5 


5000 


6 


8000 


1.00 


3 


200 


4 


1000 


5 


3000 


6 


5000 


1.005 


3 


200 


4 


1000 


5 


2000 


6 


3000 


1.01 


3 


200 


4 


1000 


5 


1000 


6 


1000 



4.2 Flux-tube profile 

We show the profiles of electric field and monopole current in Figs. 9 - 13. 
The number of configurations is N^. = 300 for all data. The investigated ranges 
are r/ro = 0.469 - 2.35 at /? = 0.98, r/ro = 0.422 - 2.53 at (3 = 0.99, 
r/ro = 0.353 - 2.12 at /3 = 1.00, r/ro = 0.305 - 1.83 at /3 = 1.005 and 
r/ro = 0.217 — 1.30 a,t P — 1.01. At glance we find that all data are clean 
enough, which allow us to identify the profile. 

For all (3 values we find a tendency that as r increases, the peak of the electric 
field at p ~ decreases and the strong peak of the monopole current profile for 
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a small p disappears. This feature can be understood as follows. For a small 
r, the mid-plane is close to the static charges so that the Coulombic electric 
field contribution is still large. In order to cancel such a strong electric field as 
much as possible, the monopole supercurrent must be high, which is indicated 
by the strong peak in Figs 9 - 13. For larger r the Coulombic field is weaker in 
the middle and consequently the peak of the monopole current is also weaker. 
It is interesting that although the rotational invariance does not hold for the 
monopole current profile, especially for small r, it is effectively restored as r 
increases. 
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Fig. 9. Profiles of the electric field (left) and of the monopole current (right) for 
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4-3 Analysis 



In this subsection, we investigate how the width of the flux-tube proflle de- 
pends on r based on the DGL analysis. In particular, we pay attention to the 
data ior r/a — 5. Corresponding flux-tube lengths are r/ro = 1.09 ~ 2.35. The 
number of conflguration used for this analysis is A^conf = 400, 500, 700, 900 and 
1200 for P = 0.98, 0.99, 1.00, 1.005 and 1.01, respectively. We have increased 
iVconf for larger (3 values to compensate the enhancement of statistical errors 
due to the smaller lattice spacings when the reference scale is introduced. 

Before the analysis, let us put all proflles into one flgure by introducing the 
Sommer scale tq. In Figs. 14 and 15, we show the proflles of the electric fleld 
and of the monopole current, respectively. We flnd that the tail of the electric 
fleld proflle (p/rg > 0.4) seems to fall into one curve, indicating its r indepen- 
dence. Moreover, although the monopole current proflle shows squeezed shape 
for small r, it becomes wider as increasing r and seems to converge again 
into one curve. Therefore, in the following DGL analysis we particularly pay 
attention to the tail behavior of the proflles. 

We briefly describe the DGL theory and the classical flux-tube solution used 
for this analysis. The DGL Lagrangian density is given by 

JCdgl = -\{d,B, - d.B^ - e*^,,f + \{d,+igBM^-\{\x?-v^f , (21) 



where and x = 0c*'' (0; f] &'^) are the dual gauge fleld and the monopole 
fleld. S^j, describes an external electric Dirac string sheet, which is bounded 
by the electric current S^S^j, = j^. This singularity is responsible for the lo- 
cation and length of the flux tube, which determines the singular part of the 
dual gauge fleld [24] . The masses of the dual gauge boson and monopole can 
be expressed as uib — \f2gv and = 2\/Xv, respectively. The inverses of 
these masses are corresponding to the penetration depth and coherence length, 
which characterize the width of the flux tube. Note that the electric cou- 
pling e and the dual gauge coupling g satisfy the Dirac quantization condition 
eg — 27r. We consider a translational invariant flux tube along the z axis by pa- 
rameterizing the system with cylindrical coordinate {p,Lp,z). As mentioned, 
due to the Dirac string E^,^ in the dual fleld strength, the dual gauge fleld 
consists of a regular and a singular parts as B^ = B^^^ + Yoi the given 

system, each part is reduced to B'^'^^ = B{p)/p and B^^'^^ = —l/{gp) e^. 
The fleld equations for B{p) and (f){p) are then given by 
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— + -— - (p - 2X(p{(p -V ) = 0. 23 

dp^ pdp \ p J 

The second term of the first equation is identified as the azimuthal monopole 
current k — k{p)e^ with k{p) — —2g(j)'^{gB — l)/p. One can solve the field 
equations analytically at large radius p where ~ f . In such region, the 
second equation provides the boundary condition of -B as -B ^ 1/g. By writing 
p = uibP and B = l/g — pK{p), the first equation can be rewritten as 



d^K IdK ( l\ 

The solution is the first-order modified Bessel function K — Ki{p). For large 
p it behaves as Ki{p) ~ y^^"^- Using this, one finds the solution for the 
electric field and the monopole current as 



E.{p) = -f- = mlMp) , (25) 
p dp 

k{p) = m%K^{p) . (26) 

We use these functions to find niB- However, we must keep in mind that 
this solution is applicable only the region where the system is translational 
invariant and the monopole field has a vacuum expectation value ~ v. 

We have employed a similar fitting procedure as used in the potential fit. 
Before carrying out the fitting, the {P*PO), where O denotes a local operator, 
and {P*P) have been averaged in bins over intervals between 20 and 40 time 
units (representing 2000 and 4000 iterations) depending on /3 values to reduce 
the autocorrelation. Then, the cylindrical profile composed only from the on- 
axis data has been computed. The mean of the dual gauge boson mass has 
been determined from the minimum of defined with the diagonal part of the 
covariance matrix, since the correlation among different p's was not significant, 
which may be due to the cylindrical averaging. The error has been estimated 
from the distribution of the jackknife samples of the fitting parameters. The 
fit range has been fixed so that the mean and/or the order of is stable 
against the change of the fit range. 

In Table 7, we summarize the result. The minimum radii which satisfy the 
above condition are found to be Pmin/o = 3 for f3 = 0.98 and 0.99, Pmin/o. = 4 
for /3 = 1.00 and 1.005, and Pmin/o = 5 for /5 =1.01. In Fig. 16, we plot m^ro 
as a function of r/ro. We find that while the mass extracted from the fit for 
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small r is larger than that from fit, it approaches Ez^s result with increas- 
ing r. Basically, if the ansatz for the DGL flux-tube solution (translational 
invariance along z axis and 4> ^ v) is valid, the masses extracted from E^ and 
those from fits should coincide with each other. In this sense we should 
take only the result for r/ro > 1.77 seriously. In fact, for short distances k^p 
cannot be translational invariant, because it is responsible for the solenoidal 
electric field inside a fiux tube, which cancels the Coulombic field at large p 
region (see. Appendix A). We find that the mass is stable around m^ro ~ 4.0. 
If the width of the flux tube diverges according to the prediction of the effec- 
tive bosonic string theory, the dual gauge boson mass, identifled within this 
analysis, should go to zero (penetration depth goes to inflnity). However, we 
have observed an almost constant behavior of the width of the flux tube in 
this range. Whether a logarithmic growth of the width is hidden in our data, 
especially for the data at r/ro = 2.35, is however difficult to say. 

Finally, we would like to discuss further the detailed flt to investigate all the 
three parameters in the DGL theory. We have performed the full range proflle 
flt (whole p region but only the on-axis data) using the flnite length flux-tube 
solution obtained numerically within the 3D lattice discretized DGL theory as 
in Ref . [25] . Here, both the electric field and the monopole current profiles have 
been fitted simultaneously, where the diagonal part of the covariance matrix 
has been taken into account to define the x^- The procedure to estimate the 
error has been the same as above. We have found, however, a tendency that 
this method cannot rehably be apphcable for P — 0.98 — 1.005, where the fine 
structure of the monopole current around the peak is not clear due to the large 
lattice spacing. In these cases, we could not identify the minimum of within 
the three parameter space, especially along the axis of the monopole mass. 
Although the dual gauge coupling has been almost a constant Pg ~ 0.06, the 
dual gauge boson mass and the monopole mass have been correlated to each 
other. The problem is that many sets of mass parameters have reproduced the 
profile well and due to this we could not find an unique set of parameters. In 
fact, for instance, as seen in Fig. 3 of Ref. [26], the monopole mass is responsible 
for the shape of the monopole current profile in small p region, while the dual 
gauge boson mass that in large p region. Only the monopole current profile at 
p = 1.01 (see. Fig. 13 or Fig. 15 at r/ro = 1-09) has the peak at p/a = 2 > 1. 
In this case we could find a minimum of and the parameters were found to 
be Pg = 0.061(2), rriBa = 0.59(2) and m^a = 0.59(6). Taking into account the 
fact that the lattice spacing is still large at P = 1.01, we may say that these 
masses are consistent with the glueball masses in the axial- vector and scalar 
channels given in Ref. [19]. 
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Fig. 14. Profiles of the electric field from different r. 
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Table 7 

The dual gauge boson mass extracted from the fit. 



/3 




nt range {p/a) 


0.98 (E,) 


1.85(2) 


3-6 


0.98 {k^) 


1.81(2) 


3-6 


0.99 (E,) 


1.75(2) 


3-6 


0.99 (k^) 


1.75(2) 


3-6 


1.00 (E,) 


1.45(2) 


4-6 


1.00 (A^v) 


1.46(2) 


4-6 


1.005 (E^) 


1.29(2) 


4-6 


1.005 (fc^) 


1.35(6) 


4-6 


1.01 


0.990(11) 


5-6 


1.01 {k^) 


1.22(2) 


5-6 
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Fig. 16. Dual gauge boson mass extracted from the fit as a function of r. 
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5 Summciry 

We have successfully simulated the Wilson gauge action of 4D compact U(l) 
lattice gauge theory on a 16^ lattice at f3 = 0.98, 0.99. 1.00, 1.005 and 1.01 
(confinement phase) by using the multi-level algorithm. 

First, we have measured the static potential and force between two static 
charges from the Polyakov loop correlation function (PLCF) up to the distance 
fma.x/f^o = 2.82. It was possible to identify the PLCF which take values from 
10~^ to 10~^^ within 10 % error. We have analyzed the potential and force 
by fitting with several ansatze and by comparing with Eq. (1) up to 0(l/r^) 
corrections, F — dV/dr — a — ^/r^. We have found that the potential ansatz 
including 7/r describes the data well and the force is consistently described by 
such a function, which have supported the universality of the 7/r correction 
to the static potential. Remarkably, we have also found that the function 
F = a — fits the force data down to relatively short distances r/ro ~ 0.3, 
which is in common with non-Abelian gauge theories. 

Secondly, we have measured the U(l) fiux-tube profile (the electric field and 
the monopole current) via the PLCF for the distances r/rg = 0.217 — 2.35. We 
have investigated the width of the flux-tube proflle as a function of r based 
on the dual Ginzburg-Landau (DGL) analysis; we have fitted the tail (large 
p region) of the U(l) flux-tube profile by the classical flux-tube solution of 
the DGL theory and have extracted the dual gauge boson mass m^, whose 
inverse characterizes the width of the fiux tube. We have found that the mass 
is almost constant for larger r as msro ~ 4.0, which indicates that the width 
remains constant in the range r/ro = 1.77 — 2.35. If we suppose that the 
width grows logarithmically as predicted by the string model, the question 
may be whether we can identify such a behavior within the limited range of r. 
This would require a further detailed investigation especially for the range 
r/ro > 1-77. 
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A composite structure of the U(l) flux tube 



In this appendix, we show the numerical evidence of the composite structure 
of the U(l) flux tube. This is possible by applying the Hodge decomposition 
to the external source to identify its monopole and photon related parts and 
by measuring the corresponding profiles. Unfortunately, since the Hodge de- 
composition requires the lattice Coulomb propagator (see, Eq. (A. 2)), which 
spoils the locality of the operator, we cannot use the LW algorithm here. Thus, 
we measure the flux-tube profile induced by the Wilson loop with a small size 
as an external source instead of a PLCF to get a clear signal. How to decom- 
pose the Wilson loop into the monopole and photon parts is given below. We 
explain this by using differential form notation. 

U(l) link variables 0{Ci) can be decomposed into the electric-photon 9^^{Ci) 
and magnetic-monopole 9"^°{Ci) parts in terms of 9{C2) and n{C2) as 

e = A-^A^ = A-\d5 + 5d)e = A-^dSe + A-^5e + 27rA-Mn , (A.l) 

where we may define 

= A-^Se , 9""° = 2TrA-^Sn . (A.2) 



To get the last equality of Eq. (A.l) we have used the relation d9 = 6 + 2?™. 
Note that the monopole part depends on n while the photon part does not. 
The Wilson loop is then expressed as 

Wa = exp[i(^, j)] = exp[i(^f ^ j)] ■ exp[i(^f ^ j)] = Wph ■ Wmo , (A.3) 



where the first term of the final expression in Eq. (A.l) does not contribute 
to the Wilson loop, since we have {A~^d66, j) = {A~^66,6j) = due to the 
conserved electric current Sj = 0. In this sense, Eq. (A.3) does not depend on 
the choice of the gauge. It is known that the static potential extracted from 
Wpfi and Wmo show a Coulombic and a linearly rising behaviors [27] . 

In order to obtain the fiux-tube profile, we measure the correlation function 



where (9 is a parity-odd local operator as Eqs. (15) and (16). Now, we may 
write O — Oph + Omo- We find that if relations 

{W)^ « {Wph)o{WMo)o (A.5) 



26 



and 







(A.6) 



are satisfied (where X and Y are arbitrary operators but contain only the 
photon or monopole part), Eq. (A. 4) can further be evaluated as 



This means that the sum of the profiles from the photon and monopole Wilson 
loops provides the full U(l) profile. The validity of the assumptions, Eqs. (A. 5) 
and (A.6), will be checked in the data. 

In Fig. A.l, we show the result with the 3x3 Wilson loop aX f3 = 0.99. Here we 
have used iVconf = 1000 configurations. These profiles are measured just on the 
mid-point of the Wilson loop. In Fig. A. 2, we also show the same electric field 
profiles as in Fig. A.l, focussing on the region around ~ 0. The important 
findings here are the following. The full U(l) profile is given by the sum of 
the photon and monopole parts, which means that Eq. (A. 7) is satisfied. The 
monopole part of E'^ becomes negative beyond a certain critical radius Pc 
(in this case pc/a ~ 1.5), indicating the appearance of the solenoidal electric 
field, which cancels the Coulomb electric field from the photon Wilson loop 
at p > Pc- The corresponding schematic picture is shown in Fig. A. 3. There 
is no correlation between the photon Wilson loop and monopole current and 
only the monopole part is responsible for the monopole current profile. This is 
consistent with the fact that the solenoidal electric field is from the monopole 
Wilson loop (dual Ampere law). Hence, we can conclude that the U(l) flux 
tube has the same composite structure as the classical flux-tube solution of 
the DGL theory [26]. This result further supports the dual superconducting 
confinement mechanism (dual Meissner efi'ect) in U(l) LGT. The behavior of 
the profiles from Wph and Wmo is completely consistent with that of the static 
potential [27]. 

We find that the peak of the electric field is -^^(O) = 0.57, which is larger than 
that in Fig. 10, ^^^(0) = 0.38, obtained by using the PLCF with r/a = 3. This 
difference exhibits the effect of the finite temporal extension t of the Wilson 
loop; the spatial part of the Wilson loop provides an additional Coulombic 
electric field. In fact, as increasing t, we can observe that the profile approaches 
to the result from the PLCF, although it becomes difficult to find a clear 
signal to check such a behavior without the Hodge decomposition method of 
the Wilson loop as shown in Figs. A. 4 and A. 5. 
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Fig. A.l. Profiles of the electric field (left) and the monopole current (right) from 
correlators with U(l) Wilson loop (3 x 3) and its photon and monopole parts. 
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Fig. A. 2. The same plots as in Fig. A.l with the axis rescaled. The profile directly 
from the full U(l) Wilson loop is omitted. 




Fig. A. 3. The composite structure of the U(l) flux tube 
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Fig. A. 4. The same plot as in Fig. A.l with W(3,5). The profiles from the full U(l) 
Wilson loop are omitted since they are too noisy. 
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Fig. A. 5. The same plot as in Fig. A.l with W(3,7). The profiles from the full U(l) 
Wilson loop are omitted. 
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